EPR entanglement strategies in two-well BEC 
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Criteria suitable for measuring entanglement between two different potential wells in a Bose- 
Einstein condensation (BEC) are evaluated. We show how to generate the required entanglement, 
utilizing either an adiabatic two-mode or dynamic four-mode interaction strategy, with techniques 
that take advantage of s-wave scattering interactions to provide the nonlinear coupling. The dynamic 
entanglement method results in an entanglement signature with spatially separated detectors, as in 
the Einstein-Podolsky-Rosen (EPR) paradox. 



One of the most important questions in modern 
physics is the problem of macroscopic spatial entan- 
glement, which directly impinges on the nature of re- 
ality. Here we analyse how rapid advances in Bose- 
Einstein condensation (BEC) in ultra-cold atoms can 
help to resolve this issue. Recently, the observa- 
tion of spin-squeezing has shown that measurement 
beyond the standard quantum limit is achievable 

Spin squeezing is known to demonstrate entan- 
glement between atoms Q, but not which subsystems 
have been entangled. An important step forward be- 
yond this would be to realise quantum entanglement 
in the Einstein-Podolsky-Rosen (EPR) sense; that is, 
having two spatially separated condensates entangled 
with each other This is an important milestone 
towards future experiments involving entanglement of 
macroscopic mass distributions, thereby demonstrating 
quantum Schroedinger cat type superpositions of dis- 
tinct mass distributions. 

In this Letter, we analyse some achievable entangled 
quantum states using a two- well BEC, and the measur- 
able criteria that can be used to signify entanglement. 
The types of quantum state considered include number 
anti-correlated states prepared using adiabatic passage, 
as well as dynamically prepared spin-squeezed states. In 
particular, we focus on spin-entanglement, as a particu- 
larly useful route for achieving measurable EPR entan- 
glement, without requiring atomic local oscillators. We 
note that spin orientation is easily coupled to magnetic 
forces to allow superpositions of different mass distribu- 
tions, once spin entanglement is present. We consider 
different types of spin entanglement criteria, and ana- 
lyze which quantum states these are sensitive to. 

We show that existing experimental techniques ap- 
pear capable of generating spatial entanglement, with 
relatively minor changes. There are several possible 
routes available. Our most significant conclusion is that 
the criterion used to measure entanglement must be cho- 
sen carefully. Not all measures of entanglement are 
equivalent, and there is an important question as to 



what one regards as the fundamental subsystems, ie, 
particles or modes. The appropriate choice of measure 
depends on the entangled state, how it is prepared, and 
what type of detection is technologically feasible. To 
demonstrate and analyse this need to adapt the crite- 
rion to the state, we choose here to analyse two and four 
mode models of a BEC, indicated schematically in Fig. 
[TJ where a±, a 2 are operators for two internal states at 
A and b\, bi are operators for two internal states at B. 





Figure 1. (a) Two internal modes a, b with spatial entangle- 
ment; (b) two pairs of modes ai, ai and bi, 62 are entangled. 



In the limit of tight confinement and small numbers 
of atoms, this type of system can be treated using a 
simple coupled mode effective Hamiltonian, of form: 
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Here n is the inter-well tunneling rate between wells, 
while gij is the intra-well interaction matrix between 
the different spin components. 

Adiabatic preparation: We first consider two-mode 
states having a single spin orientation, with number 
correlations established using adiabatic passage in the 
ground state. This makes them practical toprepare fol- 
lowing earlier experimental approaches P, [y| , as shown 
in Fig. [33a). A recent multi-mode analysis shows that 
effects of other spatial modes may be relatively small @] ■ 
In a two-mode analysis, we assume that a\ and 61 have 
been prepared in the many-body ground state of Eq ([lj 
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with a fixed number of atoms N, while the second pair 
of spin states ai and 02 remain in the vacuum state, so 
that we can write a = a\ and b = b\. In these cases there 
is only one nuclear spin orientation, and there is exist- 
ing experimental data on phase coherence and number 
correlations [3, @|, with lOdB relative number squeez- 
ing being maximally indicated. A number of previous 
analyses have used entropic measures specific to pure 
states to study entanglement. These signatures cannot 
be readily measured, and are not applicable to realistic 
mixed states that are typically created in the laboratory. 

However, one generally demonstrate spatial entangle- 
ment between the two wells a and b using the non- 
Hermitian operator product criterion of Hillery and 
Zubairy (HZ) [8J. This is also related to a recently de- 
veloped continuous- variable Bell inequality criterion Q . 
A sufficient entanglement criterion between A and B is 
the operator product measure: 



for effects of finite temperatures by assuming a canon- 
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Interwell spin operators have already been measured 
in this environment. These are defined as: 



J AB 



itfe - aSt 



/a. f± - fX +„■ fY 
/ z > J AB — J AB ICI " J 



AB' 



/(2i) 
Nab = N A 



J x 

J AB 
J AB 



Nu = 



fatS + aSt 
(at a - StS 
a) a + StS. 

In spin language, the HZ criterion shows that spatial 
entanglement is proved for any state when 
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This has similarities to the spin squeezing criterion [10( 
which has now been measured experimentally [HTD]. 
However, a crucial difference is that the spatial en- 
tanglement criterion ([3} involves an increased relative 
number fluctuation, rather than the reduced relative 
number fluctuations found with the spin-squeezing cri- 
terion. Theoretically, we find that two-well entangle- 
ment exists in the ground state with the HZ criterion, 
although suppressed for increasingly strong repulsive in- 
teractions. This behaviour is also known from previ- 
ous studies using an entropic e(p) entanglement mea- 
sure [B|, [nj . The strongest theoretical entropic entan- 
glement is found when all atom numbers are equally 
represented in the superposition. We find that the clos- 
est state to this 'super-entangled' limit is obtained at 
a critical value of Ng/n ~ —2. This attractive inter- 
action regime (as found in 41 K and 7 Li isotopes) gives 
rise to a maximal spread in the distribution of num- 
bers in each well. Maximum entanglement results for 
this model have also been found for entropic en- 
tanglement measures. In our calculations, we account 
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coupling of hn/ks = 50nK. Our result for the Hillery- 
Zubairy operator product signature is graphed below. 
This shows that two-well spatial entanglement is max- 
imized for an attractive inter-atomic coupling, and the 
effect is relatively robust against thermal excitations: 
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Figure 2. Adiabatic entanglement with interactions in a two- 
well potential. Dashed and dotted lines: HZ entanglement 
signature (E H z < 1) at T = OK, 50nK, 80nK - lowest 
lines at lowest temperature; solid line: entropic entangle- 



ment (E, 



entropic 



= 1 - e(p)/e max < 1) at T = OK. 



Dynamic preparation: To proceed further, EPR 
entanglement as we define it requires using measure- 
ments Oa and Ob that are individually defined either 
at well A or well B. Thus, entanglement is shown by 
performing a set of simultaneous measurements on the 
spatially separated systems: typically by measuring cor- 
relations < OaOb > or P(O a ,Ob)- This is necessary 
to justify Einstein's "no action at a distance" assump- 
tion, that making one measurement at A cannot affect 
the outcome of another measurement at B. One could 
achieve EPR entanglement with this criterion by mak- 
ing quadrature amplitude measurements. That is, by 
expanding a = X a + iP a and b = Xf, + iP^, where X a 
etc are "quadrature amplitudes", so that the moment 
< aw > is measured as four separate real correlations. 
Proposed methods for measuring entanglement in BEC 
experiments include time- reversed dynamics (l2^ . and 
interference with side-modes of a BEC moving through 
an optical lattice [l3| . This shows that, in principle, 
such a quadrature-based entanglement measurement is 
not impossible. However, while feasible optically, this 
type of measurement is nontrivial with ultra-cold atoms 
owing to interaction induced phase fluctuations, and we 
propose a different strategy. 

To get good EPR measurements we consider in- 
stead the intra- well "spins" J x , J Y , and J z at site 
A and B. This means having at least four modes in 
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total. To prove EPR entanglement using these mea- 
surements, one can define the spin measurements at A 

to be in terms of a\ and 0,2'- J x = (ai 02 + 02 c?i ) /2, 



alM) /2, 



aVaiJ /(2i), Jf 

#2 "2! a l so define raising and lowering 
operators as: = J x ±iJ A , and similar definition for 
site B. 

These are measurable locally using Rabi rotations and 
number measurements, without local oscillators being 
required. The spin orientation measured at each site 
can be selected independently to optimise the criterion 
for the state used. One can then show EPR entangle- 
ment via spin measurements using the spin version of 
the Heisenberg-product entanglement criterion [l4| 
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or the sum criterion by Duan et al and Simon 
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with general sum and difference spins J AB 
and J e = cos(6)J z + sin(8)J x . Here the conjugate 
Schwinger spin operators J e and J e+71 / 2 obey the un- 
certainty relation A 2 J 9 A 2 J e +*/ 2 > \\(J Y }\. 

In order to obtain ultra-cold atomic systems with 
four- mode entanglement, we consider a dynamical ap- 
proach to EPR entanglement which utilizes phase as 
well as number correlations. This requires the BEC's to 
evolve in time, in a similar way to successful EPR exper- 
iments in optical fibres [l5 , 13, 13] ■ This is very different 
to the previous scheme, as the atom-atom interaction 
appears explicitly as part of the time-evolution. The 
best entanglement is obtained when the interaction be- 
tween atoms of different spin is different to the interac- 
tion between the atoms of the same spin. In Rubidium, 
this either requires using a Feshbach resonance to break 
the symmetry, or else separating the two spin compo- 
nents spatially as in the successful fibre experiments [18( 
or in spin-squeezing atom-chip experiments At a 
Feshbach resonance, for alkali metals like Rubidium-87, 
the interactions between the different spin orientations 
can be reduced compared to the self-interactions. This 
allows an avenue for this type of entanglement with both 
the spin orientations remaining in situ in the same trap 
potential. 

To start with, we consider the conditions required 
to obtain the best squeezing of Schwinger spin op- 
erators by optimizing the phase choice 9: tg{20) = 
2{J Z , J X )/(A 2 J Z -A 2 J x ). Entanglement can be gen- 
erated by the interference of two squeezed states on a 



50 : 50 beam-splitter with a relative optical phase of 
ip = tt/2. This has also been achieved in optical experi- 
ments [l5|], although not yet in BEC. 
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Figure 3. (a) Squeezing of Schwinger spin operators 
S dB - S+ = 10/ogio [A 2 (Jl - J%)/no] (solid), S- = 

10/ogio [A 2 (J4 +7r/2 + Js +7r/2 )] /no (dashed), and n = 

\{\{ j a)\ + \(Jb)\) is shot noise (dotted), (b) Entangle- 
ment (Ep ro duct) based on the criterion Q by the solid 
curve and E 3um in sum criterion ((5| by the dashed curve. 
Here the parameters correspond to Rb atoms at magnetic 
field B — 9.131G, with scattering lengths an = 100. 4ao, 
a-22 = 95.5ao, and avz = 8O.800. ao = 53pm. The cou- 
pling constant gij oc 2coxflij- The number of Rb atoms is 
N A = 200. t = guN A t. 

Here we again take a four modes approach, explicitly 
assumingai, b\ and 02, 62 that are initially in coherent 
states for simplicity, i.e., assuming we have coherence 
between the wells. For simplicity, we suppose that the 
initial state is then prepared in an overall four-mode 
coherent state using a Rabi rotation. It is also possible 
to choose a constrained total particle number, but we 
have used the simplest model of coherence between the 
wells: 
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Next, we assume that the inter- well potential is in- 
creased so that each well evolves independently. Finally, 
we decrease the inter-well potential for a short time, 
so that it acts as a controllable, non-adiabatic beam- 
splitter to allow interference between the wells, fol- 
lowed by independent spin measurements in each well. 



4 



For dynamics, we assume a simple two-spin evolution 
per well, which is exactly soluble. We can treat this us- 
ing either Schroedinger or Heisenberg equations of mo- 
tion. In the Heisenberg case, since the number of parti- 
cles is conserved in each mode, this has the solution: 



hi (i) = exp [-i 9ijNjt a, (0) 



(7) 



where the couplings gij are obtained from the known Rb 
scattering lengths at a Feshbach resonance. 
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Figure 4. Same as Fig. [3] but assuming NO cross-couplings, 
i.e., g 12 = 0. 

After dynamical evolution from an initial coherent 
state, we find spin-squeezing in each well, prior to using 
the beam-splitter as shown in Fig. (3Ja). 

After using the beam-splitter, entanglement can be 
detected in principle as E < 1, as shown in Fig. \5j[b). 
Note, Fig. |4] shows that assuming NO cross-couplings, 
he., gi2 = gives much better results than using the 
cross-couplings obtained in a Rb Feshbach resonance. 
As discussed in the adiabatic approach, this would re- 
quire spatially separated wells for the different spin- 
orientations, or a different type of Feshbach resonance, 
in order to eliminate cross-couplings. 

In summary, we have shown two feasible techniques 
for measuring EPR-type spatial BEC entanglement, 
using currently available double-well BEC approaches 
combined with available atomic detection methods. The 
simplest method employs an attractive ground-state 
adiabatic method, with a single spin orientation. This 
requires an essentially nonlocal detection strategy, in 



which the two BEC's are expanded and interfere with 
each other. To obtain a spatially separated EPR en- 
tanglement strategy, appropriate for investigating ques- 
tions of local realism, we propose a four-mode, dynami- 
cal strategy that employs two distinct spin orientations 
in each spatial well. 
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